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Abstract 

We examine exhaustively the behavior of avalanches in critical height sandpile 
models based in two- and three-dimensional lattices of various topologies. We 
get that for two-dimensional lattices the spatial and temporal distributions 
characterizing bulk avalanches do not depend on the lattice topology. For 
the three-dimensional case, we detect a small dependence of the topology for 
the temporal distribution, while the spatial ones are independent. The two- 
dimensional lattices studied are: the plane (i?^), the cylinder {S^ xi?), and the 
Mobius-strip (M); and the three-dimensional are: , x R'^ , x x R, 
M X R, X R, K X R, and RP x R, where K and RP are respectively the 
Klein bottle and the real projective plane. 
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The notion of self-organized criticality introduce by Bak, Tang, and Wiesenfeld [|l| has 
been intensively used in the study of a huge class of phenomena in recent years. They 
showed that certain extended dissipative dynamical systems naturally evolve into a critical 
and stationary state, with no intrinsic length or time scale and independent of initial condi- 
tions, through a self-organization process. Their models are very simple simulations of real 
avalanches in a sandpile a cellular automaton in which the basic variable h{r) is an 

integer describing the height or the slope of a "sandpile" at the point r that belongs to an 
n-dimensional lattice. When /i(r) exceeds a critical value k there is an avalanche and h{r) is 
updated according to some automaton rules. This process is repeated until h{r) is less than 
k for all points r in the lattice. Many numerical and analytical studies were based in these 
models and much information is now available The scale-invariant behavior automat- 
ically generated by the very simple dynamics of these models is the main motivation to the 
intensive activity of last years in the subject. It is widely believed that the understanding 
of such toy models can shed some light in the dynamics of real systems. For instance, the 
original models of Bak, Tang, and Wiesenfeld were introduced as an intent to explain the 
ubiquitous in nature "1/f" noise [0. 

In is argued that near the self-organized critical state the system has an universal 
behavior characterized by some critical exponents. Evidences for universality in the self- 
organized critical behavior of avalanches appeared first in the work of KadanofF, Nagel, Wu, 
and Zhou [^. They studied one- and two-dimensional lattices of various sizes and their 
result pointed toward the existence of certain universality classes. For the two-dimensional 
case they also considered models with various internal symmetries, and showed that finite- 
size effects in these models are rather subtle. Scalling properties and finite-size effects for 
various types of boundary conditions were also considered P,p!0|. Robustness is another 
crucial requirement necessary to the success of these concepts in explaining real systems. 
The existence of self-organized critical states can not depend tightly of especial features of 
the model. In []T| it is shown that self-organized critical states in sandpile models are stable 
under some sort of quenched randomness. The behavior of avalanches in a random lattice 
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was considered in |]TT|, and again one gets evidences for the existence of a self-organized 
critical state, although with other critical exponents. Self-organized criticality was also 
verified under lattice cyclicities ||T2[. 

The dynamical process of these models is local in essence; when a toppling occurs in a 
given point only its nearest neighbor are affected. In spite of this, global properties of the 
lattice can alter qualitatively the time evolution of the system, as one can see, for instance, 
by comparing the number of "sand grains" that drop off the edges for lattices of different 
topologies. Indeed, one does not have a priori motivations to expect that sandpile models 
based on lattices of distinct topologies evolve in the same way for large times. In spite of 
the intense activity in this subject, the behavior of avalanches in lattices with non-trivial 
topology has not been considered yet. 

In this work we discuss the influence of global properties of the lattice on the critical 
exponents of bulk avalanches. We consider two- and three-dimensional lattices of distinct 
topologies and our result is rather intriguing: for all cases we verify the existence of a self- 
organized critical state, and we demonstrate, within our accuracy, that the critical exponents 
for the spatial distributions characterizing avalanches depend only on the dimension of the 
lattice and do not on its particular topology. As to the temporal distributions, we detect 
a small dependence of the critical exponents on the topology for three-dimensional lattices, 
while for two-dimensional ones they are independent. Our result strongly corroborates the 
robustness conjectures since we are dealing with lattices of very different global structure; 
for instance, some of them are orientables and others are not. In our models, avalanches 
are characterized by cluster sizes (s), number of topplings (m), and by fluctuation lifetimes 
(t). We have examined the following two-dimensional lattices: the plane R^, the cylinder 

X R, and the Mobius-strip M; and the following three-dimensional ones: R^, S'^ x i?^, 

X X R, M xR, S'^ xR, K xR, and RP x R, where K and RP are respectively the Klein 
bottle and the real projective plane [|13|. We simulate the distinct topologies by imposing 
certain boundary conditions on the standard lattice. Our result can be presented also as the 
fact that such boundary conditions do not alter the behavior of avalanches. Incidentally, a 



recent work [|T0| also appreciates the possibility that the presence of closed borders would 
not alter the critical behavior of avalanches in abelian sandpile models. 

We consider critical height sandpile models. In such models, the number of sand grains 
are given by the integers h{r) and the dynamics is given by the following automaton rule, 
valid for any lattice dimension d, 

if h{r') >d-l, then h{r) = h{r) - Ar'(r), (1) 

where Ar'(r) is the Laplacian matrix: Ar'(r') = 2d, Ar'(r) = —1 if r is a nearest neighbor 
of r', and Ar'(r) = otherwise. The choice of A; = d — 1 in ([^) is not important since a 
shift in k simply implies a shift in h. The rule (||) is repeated until h < d — 1 for all points 
in the lattice. If a sand grain reaches a boundary d, it shall drop off the edge. Such "free" 
boundary condition guarantees that the system stops in a minimally stable state. In order 
to impose the free boundary conditions, we assume that for r' G 9 the Laplacian matrix is 
the same one used for the bulk, although the set of nearest neighbor of r' is reduced. The 
different lattice topologies are defined from some border identification of the usual lattice, 
see Fig. In our case, they are actually imposed in the Laplacian matrix as we will explain 
bellow. 

We begin with the two-dimensional case. The usual rectangular lattice with free bound- 
ary conditions corresponds to the topology R^. The cylindrical lattice can be defined by a 
periodic condition on the standard N x M lattice: 

/i(aei + (6 + M)e2) = /i(aei + 6e2), (2) 

where ei and 62 are local generators of the lattice, and 1 < a < and 1 < b < M are 
integers. On the other hand, the Mobius-strip topology is defined by 

h{aei + {b + M)e2) = h{{N + 1 - a)ei + be^). (3) 

In the first case, if a sand grain reaches the point (a, M + 1) it will appear in (a, 1), and 
in the second one it would appear in (A^ + 1 — a, 1). In the same way, if it reaches the 



point (a, 0) it will appear in (a, M) for the cylindrical lattice and in (A^ + 1 — a, M) for the 
Mobius-strip. Such conditions can be easily implemented by means of the Laplacian matrix. 
For the cylindrical case, for instance, we need only to define as nearest neighbor of the point 
aei + Me2, 1 < a < A^, the four points (a± l)ei + (M — l)e2 and (a±l)ei + e2. Analogously, 
the nearest neighbor of aei + 62 would be (a ± l)ei + 2e2 and (a ± l)ei + Me2. In two 
dimensions one does not have others topologically inequivalent lattices with relevance to 
our purposes. The convergence of (|l]) requires free boundary conditions and thus the lattice 
cannot be compact. The only non-compact lattices that can be constructed by imposing 
boundary conditions on the standard two-dimensional lattice are the cylindrical and the 
Mobius-strip ones, see Fig. |I|. For the three-dimensional case we are also restricted to non- 
compact lattices to guarantee the convergence of (|I]). The non-compact three-dimensional 
lattice that we examine are: R^, x i?^, M x R, x x R, S'^ x R, K x R, and RP x R. 
The first three topologies are straightforward extensions of the two-dimensional cases and 
the last ones are obtained by imposing boundary conditions as follows. For an L x A^ x M 
lattice aei + he2 + 063, we define L two-dimensional sections by fixing a and considering h 
and c as free. For each one of these sections we define the Laplacian matrix according to 
Fig. |I|, imposing the boundary conditions in an analogous way of and (|^). Obviously, 
the X R topology requires A^ = M. 

In our simulations we start with a void lattice {hijc) = 0), choose at random a point r' 
and set h{r') ^ 1. The system then evolves following (|l|) until it reaches a stable state. 
This state is then locally perturbated by adding a sand grain in an aleatory point r". If the 
stability is broken, the system will evolve and we measure the "cluster size" s, given by the 
number of sites toppling at least once, the total number of topplings m, and the "fiuctuation 
lifetime", given by the number of iteration t necessary to reach again a stable state. The 
distribution of cluster sizes, of topphng numbers, and the distribution of lifetimes weighted 
by the temporal average of the response {s/t) obey respectively the following power-laws 

D{s) oc s-^, 
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Dim) oc m 



A{t) = -D{t) oc t-\ (4) 

The Table | presents the values of the critical exponents for the various topologies. Figures |^ 
and show the distribution of cluster sizes and the weighted distribution of lifetimes for 100 x 
100 and 30 x 30 x 30 lattices respectively. The distributions are averaged for each topology 
over 400 distinct samples. The situation is unchanged for lattices of different sizes. Our 
values for the spatial critical exponents (r and Tm) are in good agreement with the original 
results for two- and three-dimensional critical height models 0,0. For all distributions, the 
curves are linear over some decades but deviate for small and for large values of s, m, and t. 
As it is well known in the literature, the first can be understood as effect of the discreteness 
of the lattice and the second as finite-size effects. Also, the distribution of lifetimes gives 
less impressive results because typically we have t <C s ~ m, reducing considerably the 
number of events. However, we do not attribute the discrepancies in the 6-exponent for 
three-dimensional lattices to such a fact, we instead attribute them to genuine dependence 
on the lattice topology. 

Both set results are compatible with the condition r = r^. Such condition was first 



noticed in the numerical analysis of ||14| and it was also considered recently in ||l0[. We 
notice also that our results agree with the theoretical predictions for r and h of Zhang's 
continuous-energy model ||^, according to which tt = 1 and 6t = 1/2 for d = 2, and 
tt = 4/3 for (i = 3, although Zhang's model are based in usual hyper-cubic lattices. 

We finish noticing that, as it was already said, the distributions characterizing border 
avalanches do depend on the lattice topology. We denote by / the number of sand grains that 
drop off the edges. By applying standard finite-size scaling analysis, we fit the /-distribution 
D{f, L) for a lattice of linear size L as 

Dif, L) = L-i'g{f/Ln, for /, L » 1, (5) 

where /3 and u are the critical indices and g is the scalling function. We could verify that 
both critical indices and the scalling function do depend on the lattice topology. 
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To summarize, we have demonstrated by numerical simulations for two- and three- 
dimensional critical height sandpile models that the spatial distributions characterizing bulk 
avalanches do not depend on the lattice topology. For the temporal distribution, only for 
the three-dimensional case there is a dependence of the topology. The theoretical prediction 
of Zhang's continuous-energy model [^, tt = 2 — 2/(i, seems to be valid also with non- 
trivial lattice topologies, in spite of Zhang's original calculations were made exclusively for 
i?"'. It would be worth to understand analytically such a coincidence. Simulations on huge 
lattices, in the line of [jl4| for instance, would be of great interest to improve the data for 
the evaluation of the 6-exponent for three-dimensional lattices. 

This work was supported by CNPq, under grant 201630/93-1. N. Armesto and G. Parente 
are acknowledged for helpful discussions. 
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TABLES 





T 




h 




1.02 


1.02 


0.49 


S^xR 


1.04 


1.02 


0.52 


MxR 


0.99 


1.00 


0.50 




1.31 


1.31 


0.88 


X R"^ 


1.30 


1.31 


0.78 


S^xS^ xR 


1.31 


1.31 


0.62 


MxR 


1.29 


1.30 


0.79 


S^xR 


1.30 


1.29 


0.73 


RPx R 


1.29 


1.31 


0.77 


K xR 


1.32 


1.32 


0.75 



TABLE I. Critical exponents for avalanches in critical height sandpile models based in two- 
and three-dimensional lattices of various topologies. All exponents have an accuracy of ±0.01. 
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FIGURES 
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FIG. 1. Possible topologies obtained by imposing boundary conditions on a two-dimensional 
lattice. The sides are identified according to the arrows type and direction. Lattices (a), (c), and 
(f) are orientables, while (b), (d), and (e) are not. Only (a) and (b) are non-compact. 
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FIG. 2. Two-dimensional lattices: (a) Distribution of cluster sizes, (b) Weighted distribution 
of fluctuation lifetimes. Results for 100 x 100 lattices, each topology is averaged over 400 samples. 
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FIG. 3. Three-dimensional lattices: (a) Distribution of cluster sizes, (b) Weighted distribution 
of fluctuation lifetimes. Results for 30 x 30 x 30 lattices, each topology is averaged over 400 samples. 
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